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^ C"| Abstract. We have investigated, in the frame work of the transport approach, different aspects 

Cl^ of the QGP created in Heavy Ion Collisions at RHIC and LHC energies. The shear viscosity rj 

JL . has been calculated by using the Green-Kubo relation at the cascade level. We have compared 

the numerical results for r\ obtained from the Green-Kubo correlator with the analytical formula 
, in both the Relaxation Time Approximation (RTA) and the Chapman-Enskog approximation 

(CE). From this comparison we show that in the range of temperature explored in a Heavy Ion 
collision the RTA underestimates the viscosity by about a factor of 2, while a good agreement 
is found between the CE approximation and Gree-Kubo relation already at first order of 
approximation. The agreement with the CE approximation supplies an analytical formula that 
allows to develop kinetic transport theory at fixed shear viscosity to entropy density ratio, r/fs. 
We show some results for the build up of anisotropic flows V2 in a transport approach at fixed 
shear viscosity to entropy density ratio, rj/s. We study the impact of a T-dependent r]/s(T) on 
the generation of the elliptic flows at both RHIC and LHC. We show that the transport approach 
provides, in a unified way, a tool able to naturally describe the V2{pr) in a wide range of pr, 
including also the description of the rise and fall and saturation of the V2(pr) observed at LHC. 
Finally, we have studied the evolution of the quark-gluon composition of the plasma created 
in ultra-Relativistic Heavy Ion Collisions (uRHIC's) employing a Boltzmann-Vlasov transport 
approach that include: the mean fields dynamics, associated to the quasi-particle model, and the 
elastic and inelastic collisions for massive quarks and gluons. Following the chemical evolution 
from an initial gluon dominated plasma we predict a quark dominance close to Tc paving the 
way to an hadronization via quark coalescence. 



1. Introduction 

The experimental results accumulated in the last decade in the ultra relativistic heavy ion 
collisions (uRHICs) before at RHIC program at BNL and more recently at the LHC program at 
CERN has shown that the azimuthal asymmetry in momentum space, namely the elliptic flow 
V2, is the largest ever seen in HIC [TJ |2]. The most common approach to study the uRHICs 
is the viscous Hydrodynamics at second order in gradient expansion according to the Israel- 
Stewart theory OHMS]. The comparison of the experimental measured vi with hydro dynamical 
calculations has suggested that in these uRHICs an almost perfect fluid with a very small 
shear viscosity to entropy density ratio rj/s has been created O [5]. Similar conclusions has 
been obtained also by kinetic transport theory [6j [8] . Both Hydro dynamical and transport 



approach have shown an agreement on the evaluation of the viscosity with 4nr]/s ~ 1 — 3 very 
close to the conjectured lower bound limit 77/s = 1/Att. In general, both these calculations 
show that the elliptic flow depends sensitively on the ratio rj/s. Hydrodynamics however has 
the fundamental problem of a limited range of validity in rj/s and in the transverse momentum 
Pt < 2GeV. At increasing px viscous hydrodynamics breaks its validity because the relative 
deviation of the equilibrium distribution function 5f/f eq increases with pt becoming large 
already at pt > 3T ~ lGeV . On the other hand the relativistic kinetic transport approach 
has the advantage to be a 3+1D approach not based on a gradient expansion in viscosity that 
is valid also for large r\j s and for out of equilibrium momentum distribution allowing a reliable 
description also of the intermediate px range. 

However usually kinetic theory is applied to the study of HIC starting from the microscopic 
details of the fields and cross sections and it is not discussed directly in terms of viscosity of the 
system. The search for the QGP properties, however, have shown that the shear viscosity and 
in particular the viscosity to entropy density ratio rj/s is a key transport coefficient that could 
be very close to the conjectured lower bound limit, r//s = 1/4tt. This has lead more recently 
to develop a transport approach at fixed 77/s [U HI [9] which allows to have a direct link to the 
viscous hydrodynamic language. First attempts in this direction have been already developed 
and applied to the study of the QGP dynamics using the simple expression for 77 in the relaxation 
time approximation. On the other hand such an approach ask for a knowledge of the correct 
relation between the shear viscosity r] and temperature, cross section, mass and density. 

In this proceeding we discuss different aspect of the matter created in these uRHICs within 
the Boltzmann-Vlasov transport theory. The kinetic theory at partonic level developed solves 
the following relativistic equation: 

p^d^f(x,p) + M(x)d^M(x)d^f(x,p)=C(x,p) (1) 

where f(x,p) is the distribution function for on-shell particles and C(x,p) is the Boltzmanndike 
collision integral. We notice that the Boltzmann-Vlasov transport theory distinguishes between 
the short range interaction associated to collisions and long range interaction associated to the 
field interaction, responsible for the change of the Equation of State (EoS) respect to that of 
a free gas. This last feature allows to unify two main ingredients that are relevant for the 
formation of collective flows: the Equation of State and the finite mean free path. Furthermore 
the field interaction in M(x) is associated to a microscopic description in terms of quasi-particles 
allowing to bridge the microscopic studies of thermodynamic to the dynamics of HIC. This could 
potentially lead to infer more information on the microscopic structure of the QGP from the 
rich experimental observable of HIC's. In these proceedings we will discuss a first issue related 
with the chemical composition of the QGP, i.e. the relative abundance of quarks and gluons. 

The paper is organized as follows. In Section 2, we discuss the method for calculating the shear 
viscosity r\ using the Green-Kubo relation and after a brief overview of Chapmann-Enskog and 
Relaxation Time Approximation we show the comparison between our results using the Green- 
Kubo relation and these scheme of approximation for isotropic and anisotropic cross section. In 
Section 3, we use the results obtained in section 2 to introduce a transport approach at fixed 
shear viscosity to entropy density rj/s to study the build up of anisotropic flows V2 and the 
effect of a temperature dependent rj/s(T) on the ^(pr)- In section 4, we study the evolution of 
the quark-gluon composition of the plasma created in uRHIC's employing a Boltzmann-Vlasov 
transport approach that include: the mean fields dynamics, associated to the quasi-particle 
model, and the elastic and inelastic collisions for massive quarks and gluons. Finally Section 5 
contains summary and conclusions. 



2. Shear viscosity from the Green-Kubo relation 

The transport coefficient like heat-conductivity, bulk and shear viscosity can be related to the 
correlation functions of the corresponding flux or tensor in thermal equilibrium. The calculation 
of these coefficients is based on the fluctuation-dissipation theorem [10]. Here we are interested 
to the shear viscosity n for which the Green-Kubo formula assumes the following expression |llj : 

1 r°° f 

V =- dt d 3 x(^(x, t)**»(p,t)) (2) 
J J o Jv 

where T is the temperature, n xy is the xy component of the shear component of the energy 
momentum tensor while (...) denotes the ensemble average. In this work we determine 
numerically the correlation function {tt xv (t)iv xy (0)) solving the ultra-relativistic Boltzmann 
transport equation. 

Our aim here is to solve numerically the full collision integral to evaluate the viscosity through 
the Gree-Kubo formula and compare it with the results of the RTA and CE approximation 
scheme. The particle dynamics is simulated via Monte Carlo methods based on the stochastic 
interpretation of transition [12\ [8]. The shear component of the energy momentum tensor is 
given by 

tt^(x, t) = r^(x, t)= J ^ ^/(x, p; t) (3) 

where we notice that at equilibrium the shear stress tensor is given by the energy-momentum 
tensor. In our calculation the particles are distributed uniformly in the box. Therefore for an 
homogeneous system the volume averaged shear tensor can be written as 

i N x y 

^ y (t) = -J2 E § j - (4) 

i=i ' 

the sum is over all the particles in the box. In the Green-Kubo formula the calculation of 
the shear viscosity is reduced to the calculation of the correlation function, for details of the 
calculation of {ir xy (t)n xy (0)) see [13} I14j. Performing the numerical calculation we obtain that 
{tt xv '(t)ir xy '(0)) is an exponential decreasing function with the time. We use this fact to fit the 
correlation function with the following expression, as done in several other works [14 ^ I15|. [T6 l ll7j. 

(7r^(t)vr^(0)) = (7r^(0)7r^(0))e-* /r (5) 

t is the so called relaxation time. Substituting Eq.([S]) into Eq.([2D the formula for the shear 
viscosity becomes: 

fl = Y{ir ev (0)x x *(0))T (6) 

This is the formula that we will use in our calculation to extract the shear viscosity. The 
relaxation time r is calculated performing a fit on the temporal range where the correlation 
function assume the exponential form. The error on the value of the viscosity comes from the 
error on the initial value of the correlator and the error on the relaxation time r extracted from 
the fit of the correlation function. 

2.1. Green-Kubo vs Chapman- Enskog and Relaxation Time Approximation 
Before starting with the comparison between our results for n using the Green-Kubo relation and 
the analytical result obtained in the relaxation time approximation and the Chapman-Enskog 
approach it is useful to show briefly these two methods. The difference between these two 



methods resides in the different way in which the collision integral is approximated. In the RTA 
it has been demonstrated that the shear viscosity assumes the following expression |18t [T9] : 



V ~ J q (27 ^ 3 T a (E a ) (7) 

where T a {E a ) is the so called relaxation time and it is related to the collision frequency of the 
particles. The relaxation time can easily be expressed in terms of the total cross-section atot'- 

T^iEa) = p {a to t Vrel) (8) 

We notice that in the RTA the interaction appears in the relaxation time only through the total 
cross-section. However on general physical argument the viscosity is expected to depend also 
on the momentum transfer that on average the collisions are able to produce. In the literature 
to take into account this fact the relaxation time is approximated by r^ 1 = (pat r v re i), i.e. 
substituting the total with the transport cross section at r - This is not really coming from the 
RTA as in Eq.s ([7]) and ([8]), however we will refer to it as RTA in the following. As we will see 
such an extension of the RTA can reasonably approximate the correct viscosity only for the case 
of isotropic cross section, where it is found that t\rta = 0.8(T/&tr) = l-2(T / atot) • 

The description of the CE approach is more complicate and we use the formalism recently 
developed in Ref. [20J. At first order of approximation [QsYce can ^ e written for the most general 
case of relativistic particles colliding with a generic differential cross section a(s,Q) as 

[rjsYcz = 10 T — (9) 
coo 

where h = K? J {z)/K2{z) with z = m/T and coo = 16 (w^p — z _1 + (Sz 2 )" 1 ^^ 2 ^ . The 
are the so-called relativistic omega integrals which are given by the following integral 

/„n z 3 poo 

< J = ^732/ dyy\y-l) 7 K {2zy)a tr {y) (10) 

where j = | + i (—1)* and y = y / i/2M while at r = J dQa(s, 0) sin 2 is the transport cross 
section. 

In the following discusion we compare our results for r] using the Green-Kubo relation and 
the analytical result obtained in the relaxation time approximation and the Chapman-Enskog 
approach. We perform this comparison for two different cases: first for the case of massless 
particles colliding with an isotropic cross-section <t(s,0) = do = const and second, for a more 
realistic case, of massless particles colliding with an anisotropic cross-section. For a more general 
discussion including also the case of massive particles see [13J. For this simple case of isotropic 
cross section and massless particle the CE approximation for the shear viscosity is given by the 
following relation: 

[ V Yce = 0.8— = 1.2— (11) 

0~tr &tot 

As we can see it is the same result obtained using the RTA with r^ 1 = {po~t r v re i). In the 
literature there exist also higher order calculation up to the most recent work in Ref. [2D] where 
the calculation was extended up to the 16 th order. They higher order approximation converge 
to the value 

[rj] 1 ^ = 0.845— = 1.267— (12) 

0~tr O~tot 




Figure 1. Shear viscosity rj for a massless system as a function of the total isotropic cross-section 
a and for different temperatures. The circles are the results from the Green-Kubo method, while 
the lines are the results obtained using the Chapman-Enskog approximation to the 16*^ order 
Eq.flJ. 




Figure 2. Left: Shear viscosity rj as a function of the Debye mass mo for three different values 
of the temperature T = 0.3, 0.4, 0.5 GeV blue, green and red respectively. The open simbols 
are the results obtained using the Green-Kubo relation. The solid, dashed and dot dashed lines 
refer to the RTA approach with r _1 = (patrVrei) respectively for T = 0.3, 0.4, 0.5 GeV . The 
dotted line is the isotropic limit when uid — ► oo. Right: The same as the left panel but the 
solid, dashed and dot dashed lines refer to the Chapman-Enskog approximation at firt order. 



where however the difference between the 1° and the 16*" order is about 6%. In Fig.([T|) we show 
the results obtained using the Green-Kubo formula by full circles compared to the prediction of 
CE at \Q th order, En. (P]l by lines. The error bars for the Green-Kubo calculation are small and 
within the symbols. As we can see we have a very good agreement with the analytical results 
in all the examined range of cross sections and temperatures with a discrepancy of about 2%. 

To study the more realistic case of angular dependent cross section,we choose typical elastic 
pQCD inspired cross section with the infrared singularity regularized by Debye thermal mass 



m D : 

da 9ira^ 1 ( m 



1 + ^1 (13) 



dt~ 2 (t - m 2 D ) 



where s, t are the Mandelstam variables. Such kind of cross sections are those typically used in 
parton cascade approaches [2H [22j [8j [231 HI E]- The total cross-section corresponding to Eq. 
(fT3j) is (Ttot = Q-Ka 2 s /(2m 2 D ) which is energy and temperature independent if a s and mjj are fixed. 
In Eq. (|13p the Debye mass m/j is a parameter that regulates the anisotropy of the scattering 
cross section. For this case of anisotropic cross section the situation is more complicate and it is 
not possible to have a simple analytical expression for the viscosity. In the left panel of Fig. ([2]) 
it is shown the shear viscosity r] as a function of the Debye mass at fixed total cross-section 
atot = 3m6 and for three different temperatures T = 0.3 GeV, T = 0.4 GeV and T = 0.5 GeV. 
On the left panel of Fig. ([2]), we compare the Green-Kubo results (symbols) with the prediction of 
the modified RTA with t" 1 = (patr v rel) a t first order (lines). It is evident that there is a strong 
disagreement between the two as soon asmc is far from the isotropic limit indicated by dotted 
lines. Therefore, we see in general that even if the total cross section is fixed, non-isotropic cross 
section strongly enhance the value of the viscosity 77 and such an enhancement is stronger respect 
to i] oc cr^ 1 , i.e. the modified RTA approximation. On the right panel of Fig. ([2]), we compare 
the Green-Kubo results (symbols) with the prediction of CE at first order (lines). In this case 
we find a very good agreement between the two, hence the CE already at first order is able to 
account for the correct value of r\ even if the cross section is so forward-backward peaked to cause 
an increase of about an order of magnitude respect to isotropic limit. The RTA approximation 
would severely underestimate the viscosity. We can also see that for mo ~ 8 — 10T the isotropic 
limit is recovered and both CE and RTA coincide but this is essentially the limit discussed above. 

This study and in particular the agreement obtained beetwen the Green-Kubo calculation 
and the CE approximation supply a way to develop a transport theory with fixed rj/s, as we 
will see in the next section. 



3. The parton cascade at fixed rj/s 

In this section we introduce a transport approach at fixed 77/s ratio. The motivation to introduce 
this approach is twofold: first because it is possible to make a direct comparison to viscous 
hydrodinamic approach and second, more generally, we have a tool to directly estimate the 
viscosity of the plasma valid in a wider range of rj/s and pt respect to hydro. To this end we do 
not calculate the cross section from a microscopic model but determine the local cross section 
<Jtot in order to have the wanted local viscosity. Here we illustrate the procedure considering 
the reduction of the transport approach in to a cascade where we neglect the field interaction 
in Eq.([T]). In other words we consider the simplest case of massless gas. 

Our approach is a 3 + 1 dimensional Montecarlo cascade [8] for on-shell partons based on the 
stochastic interpretation of the collision rate discussed in Ref. |12j . In the following discussion 
we use the pQCD inspired cross section used in the previous section Eq. (|13p . In the CE 
approximation for anisotropic cross section it is not possible to express the shear viscosity rj 
in an analytical form, because the differential cross section enters in the so called relativistic oj 
integrals Eq. (|10p . For this particular case 77 can be written in the following form: 

HcE=9{m D ,T)— (14) 

<?tot 

g(rriD,T) is a function of temperature and the thermal Deby mass. Therefore considering that 
the entropy density for a massless gas is s = p(4 — fJ>/T), \i being the fugacity, we get: 



(15) 




Figure 3. Different temperature dependent parametrizations for r]/s. The orange area take 
into account the quasi-particle model predictions for rj/s [25] . 



where otot 1S the total cross section. In our approach we solve the relativistic Boltzmann equation 
with the constraint that rj/s is fixed during the dynamics of the collisions in a way similar to 
|24j but with an exact local implementation as described in detail in [8]. We can evaluate locally 
in space and time the strength of the cross section o~tot(p,T) needed to have n/s at the wanted 
value by mean of the following formula: 

g(m D ,T) T 

- w=m Wce/s (16) 

An alternative but equivalent way to look at such procedure is that we implement a total 
cross section of the form atot = K(p,T)a p QCD > &pQCD where K takes into account the 
non perturbative effects responsible for that value of viscosity r](T). Note that this approach 
have been shown to recover the viscous hydrodynamics evolution of the bulk system [4, 8j, but 
implicitly assume that also high pt particles collide with largely nonperturbative cross section. 



3.1. Effect of temperature dependent rj/s(T) on V2 

In our calculation the initial condition are longitudinal boost invariant. The initial dN/dr] 
have been chosen in order to reproduce the final dN c h/dr]{b) at mid rapidity observed in the 
experiments at RHIC and LHC energies. The partons are initially distributed in coordinate 
space according to the Glauber model while in the momentum space at RHIC (LHC) the 
partons with px < Po = 2GeV (AGeV) are distributed according to a thermalized spectrum 
with a maximum temperature in the center of the fireball of To = ITq (To = 3Tc) and the 

-» -* 1/3 

transverse profile T(r) = To(p(r) / p(0)) while for pr > Po we take the spectrum of non- 
quenched minijets according to standard NLO-pQCD calculations. We also start our simulation 
at the time t$ = 0.6/m/c at RHIC and to = 0.3/m/c at LHC in agreement with the t^T ~ 1 
criterium. In order to study the effect of the kinetic freezeout on the generation of the elliptic 
flow we have performed different calculations: one kind of calculation with a constant Anrj/s = 1 
and 2 during all the evolution of the system (black and green dashed line of FigJ3]) the other 
(shown by red solid line in FigED with Airn/s = 1 in the QGP phase and an increasing rj/s in 
the cross over region towards the estimated value for hadronic matter Anrj/s = 8 |26j . Such 
an increase allows for a smooth realistic realization of the kinetic freeze-out. In Fig. [J] it is 
shown the elliptic flow V2{pr) at mid rapidity for 20% — 30% centrality for both RHIC Au+Au 
at Vs = 200GeF (left panel) and LHC Pb+Pb at </s = 2.76TeV (rhigh panel). As we can see 
at RHIC energies, left panel of Fig. [H the Vi is sensitive to the hadronic phase and the effect 
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Figure 4. Differential elliptic flow V2(pt) at m id rapidity for 20% — 30% collision centrality. 
On the left panel, the orange band indicate RHIC results measured by STAR and the orange 
points on the right panel are the LHC results measured by the ALICE collaboration, data taken 
by [2J. The black and green dashed lines are the calculations with Anrj/s = 1 or 2 respectively 
during all the evolution of the fireball and without the freeze out condition while the red line 
is the calculation with the inclusion of the kinetic freeze out and with 4tti]/s = 1 in the QGP 
phase respectively. 

of the freeze out is to reduce the V2 of about of 25%, from black dashed line to red solid line 
and the effect of the kinetic f.o. is quantitatively similar to have a constant Airrj/s = 2 during 
all the evolution of the fireball, green dashed line in Fig. SJ For the pr range shown we get 
a good agreement with the experimental data for a minimal viscosity rj/s ~ 1/(4tt) once the 
f.o. condition is included. This result is in agreement with the viscous hydro results, even if we 
notice that the agreement is good up to px ~ 3 GeV and non ansantz on the non-equilibrium 
effect is needed. At LHC energies, right panel of Fig. 01 the scenario is different, we have that 
the V2 is less sensitive to the increase of rj/s at low temperature in the hadronic phase. The effect 
of large rj/s in the hadronic phase is to reduce the v% by less than 5% in the low pr region, from 
black dashed line to the red solid line in right panel of Fig. [H This different behaviour of V2 
between RHIC and LHC energies can be explained looking at the life time of the fireball. In fact 
at RHIC energies the life time of the fireball is smaller than that at LHC energies, 4 — 5/m/c 
at RHIC against the about 8 — Wfm/c at LHC, in agreement with HBT results. Therefore 
at RHIC the elliptic flow has not enough time to fully develop in the QGP phase. While at 
LHC we have that the V2 develops almost completely because the fireball spend more time in 
the QGP phase. Qualitatively these results are similar to those obtained using the formula for 
rj in the RTA (see Ref.([27j)) but quantitatively the results are different because as shown in 
the previous section the viscosity estimated in the RTA in the range of temperature explored 
in a HIC differes also for a factor 2 from that obtained in the CE approximation which is the 
correct one. In general in these results we have a smaller V2 respect to the results obtained 
used the RTA approximation for rj and this is due to the fact that for fixed rj/s ratio in the CE 
approximation we estimate locally a smaller total cross section respect to the RTA case. 

As pointed out due to this large life time of the fireball at LHC and the larger initial 
temperature is interesting to study the effect of a temperature dependence in rj/s. In the 
QGP phase rj/s is expected to have a minimum of rj/s (47r) _1 close to Tc as suggested by 
1QCD calculation [2H]. While at high temperature quasi-particle models seems to suggest a 




Figure 5. Differential elliptic flow V2{pt) at mid rapidity for different collision centralities. 
On the left panel, the blue red and green bands indicate RHIC results measured by STAR and 
while points on the right panel are the LHC results measured by the ALICE collaboration with 
the same colours of the left panel for the different collision centralities, data taken by [2J. The 
solid and dashed lines are calculations with the inclusion of the kinetic freeze out and with 
irrrj/s = 1 and Airrj/s oc T 2 in the QGP phase respectively. The different colour are for the 
different collisions centralities. 



temperature dependence of the form rj/s ~ T a with a ~ 1 — 1.5 [25J. To analyze this possible 
scenarios for rj/s in the QGP phase we have considered another different with a quadratic 
dependence Airrj/s = (T/T ) 2 = (e/e ) 1/2 ( green line) where eo = 1.7GeV/ fm? is the energy 
density at the beginning of the cross over regions where the rj/s has its minimum, see blue 
dashed line in Figj3j At RHIC energies the V2 is essentially not sensitive to the dependence of 
rj/s on temperature in the QGP phase, see the dashed blu, red and green lines in the left panel of 
Fig. [5j However the effect on average is to decrease the value of V2 but at low pt < 1.5GeV the 
V2(pt) appears to be insensitive to r//s(T) while a quite mild dependence appears at higher px 
where however the transport approach tends always to overpredicted the elliptic flow observed 
experimentally. At LHC energies the build-up of V2 is more affected by the rj/s in the QGP 
phase and on average it is reduced of about a 20%. In any case still a strong temperature 
dependence in rj/s has a small effect on the generation of V2 we found that with a constant 
or at most linearly dependent on T rj/s the transport approach can describe the data at both 
RHIC and LHC at least up to pt ~ 2 — 3GeV. The fact that we can reproduce the V2(pr) for 
different centralities implies that we can reproduce also the breaking of Vije scaling as described 
by core-corona models [291 1301 l3T] . In our approach, this fact is taken into account dynamically 
by the increase of the rj/s ratio at lower energy density, see Fig. [3j 

3.2. Effect of high px partons on V2 

In the previous section has been pointed out that in our approach where we fix locally the ratio 
■q/s we have that axot = K(p,T)a p QCD therefore we have large cross section independently of 
the pt of the colliding particles. Obviously, this procedure doesn't permit to recover the pQCD 
limit for hard collision therefore we extend our previous approach allowing for a K factor that 
depends on the invariant energy of the collision K(s). We choose this function in such a way 
that at high energies K(s) — > 1 and in this way we get the the connection between the non 
pertubative interacting bulk and the asymptotic pQCD limit. For K(s) we choose the following 




Figure 6. Left: V2(pt) at mid rapidity for 20% — 30% collision centrality at RHIC. The dashed 
lines are the calculations with K = const and for Attv/s = 1 and Attv/s oc T with f.o. respectively 
for black and blue curves while the solid lines are the same but with K(s/A 2 ). Right: v%{pt) at 
mid rapidity and for 20% — 30% collision centrality at LHC with the same legend, data taken 
from [2]. 



exponential form K(s/A 2 ) = l + je~ s ^ 2 , where A is a scale parameter that fix the energy scale 
we have a transition to pQCD behaviour. While 7 plays the same role of K in the old calculations 
and it is determined again in order to keep fixed locally the rj/s. Therefore we can repeat the 
same procedure as described in the previous section but now with OTot = K(s)a p QCD- Only 
for a qualitative discussion in this section we will show the results for the elliptic flow obtained 
using the relaxation time approximation (RTA) for the shear viscosity [27]. Implementation 
with the CE is in progress. In the left panel of Fig. [6] we compare the V2(pt) for 20% — 30% 
collision centrality at RHIC energies with (solid) and without (dashed) the inclusion of an energy 
dependent K factor and for two T dependences of the 77/s. As we can see K{s/A 2 ) does not 
affect at all the V2(pt) for pr < 2GeV, in other words at RHIC energies the high px partons do 
not affect the generation of the V2 of the bulk. On the other hand, we have a reduction of the 
V2 for pr > 3GeV and in particular the v% becomes a decreasing function of px for pt > 3GeV 
in agreement with what is observed experimentally (orange band). In the right panel of Fig. 
[6] we compare in a large range the V2(pt) at LHC energies with (solid) and without (dashed) 
the inclusion of an energy dependent K factor and for two T dependence of the rj/s. We notice 
that the two sets of experimental data refer to different method of v% measurements, namely 
t>2[2] (circle) and ^[4] (square) and our theorethical results should be compared to t>2[4] because 
event-by-event fluctuations are not considered in our calculations. As we can see at LHC energies 
the V2 is sensitive to K(s/A 2 ) already at pt ~ 1.5GeV quite lower than the RHIC case, in other 
words the many high pr partons that we have at LHC energies affect the generation of the V2 
of the bulk. In general we observe that the V2(pt) becomes more sensitive to the value of the 
viscosity when we include the function K(s/A). 

We can give the following interpretation for the V2{pt)'- The raise of the V2{pt) at low pt 
is an effect of a strong interacting fluid where the particles interact mainly non perturbatively 
with large cross sections and therefore we get a description in agreement with hydrodynamics. 
For pt > 3 — 4GeV the elliptic flows starts to be a decreasing function of pr- This is the region 
where the disappearance of the non perturbative effect significantly affects the V2{pt) making 
faster and stronger (~ 20 — 25%) the fall in the elliptic flow in the range 3GeV < pr < 8GeV. 
Finally in our calculation for pt > 8GeV, where the pQCD limit is almost established, seems to 
appear the saturation of the V2 similarly to the experimental data and typical of a path-lengh 



mechanism as in jet quenching models |32j . In this region an analysis with better statistics is 
required. 

4. Chemical equilibration of the QGP 

A successful way to account for non-perturbative dynamics of the QGP is a quasi-particle 
approach, in which the interaction is encoded in the quasi-particle masses [221 1331 EH EHl 
E3 EH [38]. The model is usually completed by introducing a finite bag pressure that can 
account for further non-perturbative effects. It is well known that, in order to be able to 
describe the main features of lattice QCD thermodynamics [39], in these quasi-particle models a 
temperature-dependent mass has to be considered. This also implies that the bag constant has 
to be temperature-dependent, in order to ensure thermodynamic consistency. We notice that if 
the QGP can be described in terms of finite mass excitations this has a strong impact on the 
quasiparticle chemical ratio N q+q /N g . In fact at equilibrium one has: 

AW = d q+ - q m 2 q {T)K 2 {m q /T) 
N g d g m? g {T)K 2 {m g /T) ' 1 > 

where K 2 is the Bessel function and m qtg (T) are the T— dependent quark and gluon masses 
that can be determined by a fit [25] to recent 1QCD calculations [39]. In Figll71 we show 
by solid line the equilibrium ratio when the fit to 1QCD e(T) is done assuming m q /m 2 = 
3/2 • (A c 2 - 1)/A C (2A C + Nf) = 4/9 according to a pQCD scheme [351 ESI EZ] . 
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Figure 7. Quark to gluon ratio at equilibrium as a function of the temperature as predicted 
by QPM [25]. For the solid line the m q /m g ratio is according to pQCD while the dashed line 
indicate the massless case. 



The motivation to use a transport equation for quasiparticles with a space-time dependent 
mass resides in the success of quasi particles in describing correctly the behavior of energy density 
and pressure of the QGP as computed in the 1QCD approach. We employ a Boltzmann-Vlasov 
transport theory to simulate the partonic stage of the HIC. In the last years several codes have 
been developed based on transport theory at the cascade level [I0J [221 E21 IE] , i.e. including only 
collisions between massless partons, with quite rare exceptions [41[ l9~l l4~2]. 

In this section we present the results within a transport approach that includes the mean field 
dynamics associated to the thermal self-energies generating the finite mass m(T) in the quasi- 
particle model discussed in Ref . [23 E3 E3 EZ] • In such a picture the relativistic Boltzmann- 
Vlasov equation can be written as follows: 



[p^ + m*(x)d^m*(x)d»]f(x,p) = C[f](x,p) 



(18) 



where C(x,p) is the Boltzmann-like collision integral, main ingredient of the several cascade 
codes: 

C= flhvfv-f^Mvv^vtf&ipi +P2-P1-P2) (19) 



21'2' 

where Jj = Jj d 3 pj /(27r) 3 2Ej, fj are the particle distribution functions, while -Mf~>i denotes the 
invariant transition matrix for elastic as well as inelastic processes. The elastic processes have 
been implemented and discussed in several previous works [5DI [221 U21 E] • In this section, instead, 
we will show some results including the inelastic processes between quarks and gluons (gg qq) 
and to achieve this we have evaluated the matrix element in a pQCD LO order scheme. The tree 
diagrams contributing to the gg qq correspond to the u, t, s— channels: M. = M. s + Ait + -M u . 
For the massless case the cross sections for such processes are the textbook pQCD cross section 
for jet production in high-energy proton-proton collisions. In our case we have considered a finite 
mass for both gluons and quarks together with a dressed gluon propagator where for vanishing 
gluon mass we recover the renowed Combridge cross sections [33] . The details of the calculations 
are quite similar to the one in [H] for finite current strange quark mass. The thermodynamical 
self-consistency of the QPM requires a self-consitency between the Bag constant and the effective 
mass of the quasiparticles [25] which leads to a gap- like equation coupled to Eq. (fT8j) : 



dB f d 3 p rm{x) 

-fi{x,P) (20) 



drrii J (2tt) 3 Ei(x) 

with i = q,q,g. Eq. (|20p allows to evaluate locally the mass in Eq. (|18p also in non-equilibrium 
conditions guaranteeing the conservation of the energy-momentum tensor of the fluid. 




Figure 8. Left:Quark to gluon ratio as a function of time (normalized to To) in Au + Au at 
y/s = 200AGeV (black lines) and for Pb + Pb at \fs = 5.5ATeV (light lines). Dashed lines are 
for the massless case and the solid for the massive case. Right: Quark to gluon ratio vs pr at 
the freeze-out time. The thin solid line represents the full equilibrium raio, Eq. (|2ip : see text 
for details. 



We have checked that the numerical implementation in a stationary box is able to reproduce 
correctly both the kinetic and the chemical equilibrium. In particular we have carefully checked 
that for different temperature we get the correct abbundancy of quarks and gluons as given by 
Eq. (|17p . In the following discussion we will show some results for the more realistic case of a 
HIC. 



We have simulated Au + Au collisions at yfs = 200 AGeV and Pb + Pb at y/s = 5.5ATeV for 
— 10% centrality. The initial condition are longitudinal boost invariant with an initial quark 
to gluon ratio of N q+q /N g = 0.25 with about 75% of gluon and 25% of quarks. The initial 
conditions in the r-space are given by the standard Glauber condition while in the p-space we 
use a Boltzmann-Juttner distribution function up to a transverse momentum pr = 2 GeV while 
at larger momenta minijet distributions are implemented, as calculated by pQCD at NLO order 
|45j . At RHIC the initial maximum temperature at the centre of the fireball is To = 340 MeV 
and the initial time tq = 0.6 fm/c (corresponding also to the to -To ~ 1 criterium) as in succesfull 
hydro dynamical simulations. For Pb + Pb collisions at y/s = 5.5 ATeV we have To = 600 MeV 
and To ~ 1/To = 0.3 fm/c. In FigJ5]it is shown the time evolution of the ratio R qg = N q+q /N g 
for both the massless case (dashed lines) and the massive quasi-particle case (solid lines). As 
we can see in the massless case at RHIC the system reaches very quickly, in less than 1 fm/c, 
the chemical equilibrium given by R qg ~ 2. The difference in the time evolution of the ratio 
between RHIC and LHC is due to the temperature dependence of the strong coupling as(T) 
that we have considered in our model. In fact during the first stages of the fireball evolution 
at LHC the temperature is larger with respect to the initial stages at RHIC and this means a 
smaller as{T) at LHC. However for the final ratio we obtain similar results for RHIC and LHC 
and this behaviour is consistent with the fact that at the equilibrium this ratio does not depend 
on the temperature of the system, see dashed line in Fig.©. 

When the quasi-particle massive case is considered, as we can see the ratio R qg is an increasing 
function with the time and it reaches R qg ~ 3.8 for Au + Au. As we can see in Fig. ([8]) for 
the massive quasi-particle case (solid lines) we have a difference in the time evolution of the 
ratio N q+q /N g between RHIC and LHC and this is consistent, not only with the temperature 
dependence of as(T), as in the massless case, but also with the fact that the equilibrium value is 
strongly T dependent especially close to T c (see FigjT]) and the system is sensitive to the different 
initial temperature between RHIC and LHC. Nonetheless we find that the fireball reaches a value 
relatively close to the equilibrium at T ~ T c and eventually it is composed by about 80% of 
quark plus anti-quarks. As shown we have that the time evolution of the quark to gluon ratio 
N q+q /N g is qualitatively similar between RHIC and LHC but for LHC it is always smaller than 
that at RHIC. This difference is due to the fact that at LHC longer part of the lifetime is spent in 
a T-region where the equilibrium R qg is nearly constant. This results into a moderately smaller 
final ratio. 

In the right panel of FigJH] the px dependence of the quark to gluon ratio is shown for the 
initial distribution (dotted line) and the freeze-out distributions: massless case (dashed line) 
and massive case (solid line). Black lines are for Au + Au and light ones fore Pb + Pb. We see 
the large difference between the massless and the massive case and also that the net gluon to 
quark conversion extends up to quite large pt- As we can see we have that at very low px this 
ratio is very close to the equilibrium value for both massive with N q+q /N g ~ 6 — 7 and massless 
with N q+q /N g ~ 2 at the f.o. where the temperature is T ~ Tc- The decrease of the ratio with 
Pt can be expected considering that high— px particles can more easily elude the equilibration 
dynamics. However, in the massive case, we note a quite strong dependence below px = 2 GeV 
that has not to be interpreted as a fast detachment from the chemical equilibrium. In fact the 
Pt - dependence of the ratio can be evaluated analitically at equilibrium and it is given by 

dN/d 2 p T \ q+q = d q+q m q T e^ m T^T)/T] 
dN/d 2 p T \ g d g m 9 e 7[(»4-A,PT)/T] [ ' 

where (3q is the radial flow velocity, m T and m T are the transverse masses. In the right panel of 
Fig. ([8]) with the thin dashed line we plot the function of Eq, (f2"Tj) rescaled by a factor 0.85 
accounting for the lack of full thermalization. The strong px dependence obtained in the 
transport simulation follows very closely the equilibrium behavior at least up to px ~ 1-5 



GeV. This is a well known effect predicted by hydrodynamics and experimentally observed from 
SPS to LHC energy for hadronic spectra. In conclusion the effect of the mass on the chemical 
composition of the QGP is substatially both at high px and low pr region where we have a 
dominance of quarks also in high pt region, this result represent an important support to the 
idea of a coalescence of masssive quraks [Ml U5J H7J, 05J . 

5. Conclusions 

In the frame work of the transport approach we have investigated different aspects of the QGP 
created in Heavy Ion collisions at RHIC and LHC energies. 

We have developed a method to solve numerically the Green-Kubo formula for the shear 
viscosity for the case of a relativistic Boltzmann gas. In our study we have compared the 
Chapman-Enskog approximation and the relaxation time approximation with the result from 
the Green-Kubo correlator that in principle should provide the correct result. We have shown 
our results only for two physical case: isotropic and non-isotropic cross section for massless, 
for a more complete study see [13]. Our work shows that the relaxation time approximation 
always underestimate the shear viscosity even by more than a factor of 2 and it gives a 
satisfying prediction only in the unrealistic case of massless isotropic cross section. Instead 
in the Chapman-Enskog approximation already at first order have shown an agreement at the 
level of 2% with the numerical calculations of the Green-Kubo correlator for all the physical 
case considered. The agreement between the CE approximation at first order and the Green- 
Kubo method also supplies a relatively simple analytical expression for the viscosity that we 
have used to developed a kinetic transport theory at fixed viscosity with very good precision. 
This shows that the current used relaxation time approximation for the viscosity can lead to 
significantly underestimate the rj/s. Therefore previous works [8j H9| [27] in this direction are 
only approximately valid. 

We have investigated within a transport approach at fixed rj/s, in the CE approximation, 
the effect of a temperature dependent rj/s at RHIC and LHC energies. An important result is 
that at LHC a key observables like the elliptic flow is much less contaminated by the hadronic 
phase allowing a better study of the QGP properties. This result is qualitatively in agreement 
with the previus work in the RTA. We get for both at RHIC and LHC a good agreement with 
the data when the ratio rj/s ~ l/(4?r) and in general we observe not a large sensitivity of v% 
to the T dependence in r//s. Furthermore we have seen that at LHC the large ammount of 
particle with p? > ^GeV/c interacting nearly perturbatively cannot be neglected. The interplay 
between perturbative and non-perturbative behaviour seems to have an important effect on the 
generation of t>2 at intermediate px and it could explain the rapid raise and fall of V2{pt) m 
< pt < 8GeV/c shown in the experiments. 

Finally, we have studied the evolution of the quark-gluon composition of the plasma 
created in ultra- Relativistic Heavy Ion Collisions (uRHIC's) implied by standard quasi-particle 
approach developed to study the thermodynamics of the QGP. Tis has been realizd employing a 
Boltzmann- Vlasov transport theory that includes both elastic and inelastic collisions. This study 
shows that one could expect that the QGP created in uRHIC's, even if it is initially a Glasma 
should very quickly evolve into a plasma dominated by quark plus antiquarks close to the cross- 
over temperature T c . With the quark to gluon ratio can evolve by more than a factor 20 and 
at freeze-out is anyway almost a factor 2 larger than the one for an equilibrated massless QGP. 
The results are quite robust and devolpments of the QPM [251 EHJ ESI EZ] or inclusion of three- 
body inelastic scatterings may even make the effect larger. The result supplies a theorethical 
justification of the massive-quark coalescence hadronization models able to successfully describe 
several puzzling observations at RHIC and LHC [Ml [Ml H7J HH] . 
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